Dynamics of quantum spin liquid and spin solid phases in IPA-CuCla under field. 
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Inelastic and elastic neutron scattering is used to study spin correlations in the quasi-one dimen- 
sional quantum antiferromagnet IPA-CuCls in strong applied magnetic fields. A condensation of 
magnons and commensurate transverse long-range ordering is observe at a critical filed He ~ 9.5 T. 
The field dependencies of the energies and polarizations of all magnon branches are investigated 
both below and above the transition point. Their dispersion is measured across the entire ID Bril- 
louin zone in magnetic fields up to 14 T. The critical wave vector of magnon spectrum truncation 
[Masuda et al, Phys. Rev. Lett. 96, 047210 (2006)] is found to shift fi-om « 0.35 at H < He 
to he = 0.25 for H > H^. A drastic reduction of magnon bandwidths in the ordered phase [Garlea 
et al, Phys. Rev. Lett. 98, 167202 (2007)] is observed and studied in detail. New features of the 
spectrum, presumably related to this bandwidth collapse, are observed just above the transition 
field. 

PACS numbers: 75.10.Jm, 75.25.+Z, 75.50.Ee 



I. INTRODUCTION 



Recent years were marked by substantial progress in 
the study of gapped quantum-disordered antiferromag- 
nets (AFs) in external magnetic fields. The most promi- 
nent common feature is the destruction of the non- 
magnetic spin liquid ground state and the emergence of 
a long-range AF ordered phase at some critical field He. 
Such quantum phase transitions can often be viewed as 
a condensation of magnons. The Zeeman effect drives 
the energy of a particular gapped magnon to zero at 
the AF wave vector, whereupon they are macroscopi- 
cally incorporated in the ground state. By now, the 
phenomenon has been observed and studied experimen- 
tally in numerous prototype materials. In spin systems 
with relevant magnetic anisotropy, such the 5 = 1 Hal- 
dane spin chain compound NDMAP,ii^ the model bond- 
alternating S = 1 chain NTENP^ or the S = 3/2 



spin-dimer system Cs3Cr2Br9^i^ the transition is of an 
Ising universality class,^ and the high-field ordered phase 
is gapped. More recent work focused on materials de- 
scribed by Heisenberg or XY (axially symmetric) Hamil- 
tonians: SD-interacting S = 1/2 dimers in TlCuClaji^ 
the S ^ 1/2 quasi-2D network PHCC^, the 5 = 1/2 
square lattice bilayer BaCu2Si2 07,^^ and the S ~ 1 
material NiCl2-4SC(NH2)2-^^ Here the transition can 
be viewed as a Bose-Einstein Condensation (BEC) of 
magnons,^* despite the ongoing controversy regarding 
such a nomenclature.^- Correspondingly, the high field 
ordered phase remains gapless. 

In all the diverse cases, due to the soft-mode nature of 
the phase transition, the key to the underlying physics 
is in understanding the spin excitations on cither side of 
the phase boundary. Until recently, of the more isotropic 
materials, only in TlCuCla has the spectrum been stud- 
ied in sufficient detail. Unfortunately, recent data indi- 
cates that the transition in this compound is actually not 
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BEC, and that the high field spectrum is gapped, due 
to Dzyaloshinski-Moriya type anisotropyiS and magne- 
toelastic effects.— In a recent short paper we reported 
neutron scattering studied of magnetic excitations in the 
S = 1/2 spin ladder IPA-CuCla, where the transition is 
an almost exact realization the BEC case.'^^ In fact, any 
long-range properties, such the emergence of a gapless 
collective Goldstone mode beyond the transition point, 
appear to be fully consistent with original model of Gi- 
amarchi and TsvelikJ^ Interestingly, the transition was 
also found to dramatically affect short-range spin corre- 
lations and excitations near the magnetic zone boundary. 
This phenomenon was attributed to AF long-range order 
violating discrete translational symmetry. Such a spon- 
taneous symmetry breaking is absent in either the con- 
ventional BEC, or in magnon condensation transitions in 
NTENP, TlCuCla and many other spin gap systems. 

The present paper is a report on extensive high-field 
neutron scattering studies of IPA-CuCla . In addition to 
providing more details on the previous experiments and 
data analysis procedures, we focus on a number of new 
results. Among these are a study of polarizations of mag- 
netic excitations, observation of certain unusual features 
of the spectrum just beyond the transition point, a study 
of the effect of long-range ordering on single-magnon 
spectrum termination,^^ and a measurement of magnon 
dispersion at the highest attainable field of 14.5 T. 

IPA-CuCls crystallizes in a triclinic space group PI 
with a = 7.766 A, 6 = 9.705 A, c = 6.083 A, a = 97.62°, 
P = 101.05°, and 7 = 67.28°.20 As discussed in detail in 
Ref . [l9l. the key features of the structure are S* = 1/2 lad- 
ders of Cu^"*" ions that run along the crystallographic a 
axis, and approximately parallel to the (a, c) plane. Spin 
correlations are ferromagnetic on the ladder rung, but an- 
tiferromagnetic along the legs. The thus-formed effective 
S = 1 chain has a Haldane gap2i of A w 1.5 meV )^^i^° 
seen at the ID AF zone-center q = ( ^"^"'^ ,k,l), where n 
is integer. Due to weak interactions between ladders in 
the (a, c) plane, the gap varies between A ~ 1.2 mcV at 
k — (global dispersion minimum) and A = 1.8 meV 
at k = 0.5. Interactions between ladders along the b 
axis are negligible, due to entirely non-magnetic sheets 
of organic molecules in-between. At low temperatures 
IPA-CuCls goes through a field-induced AF ordering 
that manifests itself in a lambda specific heat anomaly 
and the appearance of non-zero uniform magnetization.-^ - 
The transition occurs at i?c,a = 10.2 T, Hc,b = 10.5 T 
and Hc.c = 10.4 T, in fields applied along the crystallo- 
graphic a*, b* and c* axes, respectively. This anisotropy 
effect is fully accounted for by the anisotropy of the g- 
tensor for Cu^+ in this system: ga = 2.05 gb = 2.22 and 
5^ = 2.11^^ 



II. EXPERIMENTAL 

In this work we summarize the results of several sepa- 
rate series of neutron scattering experiments on deuter- 



ated single-crystal IPA-CuCla samples. Measurements 
on the SPINS cold-neutron spectrometer at the NIST 
Center for Neutron Research (NCNR) were performed 
on an assembly of 20 single crystals of total mass 3g, 
mounted with the crystallographic b axis vertical, and 
the {h, 0, 1) reciprocal space plane accessible for measure- 
ments (Setup I) . The mosaic of the assembly was some- 
what irregular, with a FWHM of 4.5°. Scattered neu- 
trons were analyzed at 3.7 meV by a pyrolitic graphite 
PG(002) analyzer in horizontal monochromatic-focusing 
mode, used in combination with a BeO filter. A 120' 
radial collimator was installed between sample and an- 
alyzer. A magnetic field of up to 11.5 T was applied 
along the b axis. A dilution refrigerator maintained sam- 
ple temperature at 100 mK. In Setup II the sample was 
mounted with the c axis vertical (along the applied field) 
and the (/i, fc,0) plane open to scattering. An 80' coUi- 
mation was used in front of the sample, and only 5 out 
of 11 analyzer blades were focused on the detector. Ex- 
periments in Setup III were performed at the V2-FLEX 
3-axis spectrometer at HMI. A smaller crystal assembly 
of approximate mass 1 g was mounted with the b axis 
vertical in a 14.5 T magnet. The data were collected at 
70 mK, with no dedicated coUimation devices in the neu- 
tron beam, a PC filter after the sample, and a focusing 
PC (002) analyzer tuned to 3.7 meV. A 15 Tesla mag- 
net was used in Setup IV, where the data were taken at 
T < 100 mK on a 3g sample mounted with the c-axis 
vertical on the IN-14 cold-neutron 3-axis spectrometer 
at ILL. Sample mosaic was about 4° FWHM. Scattered 
neutrons were analyzed at 2.9 meV by a focusing PC 
(002) analyzer with a Be filter. A 60' beam collimator 
was used in front of the sample. 

In addition to these 3-axis data, the discussion be- 
low will involve time-of-flight (TOF) spectra previously 
collected using the Disk Chopper Spectrometer at NIST 
(Fig. 1 of Ref. 18). The sample in this experimental Setup 
V was mounted with the c axis vertical, parallel to the ap- 
plied field. Incident neutron energy was fixed at 6.7 meV. 
Finally, diffraction data were taken on the D23 lifting 
counter diffractometer, at ILL, using a monochromatic 
beam with A — 1.276 A (Setup VI). The measurements 
were carried out on a 3.2 x 2.3 x 8.5 mm'^ deuterated 
single crystal specimen. For this experiment we used a 
12 T vertical field superconducting magnet with a dilu- 
tion insert. 



III. RESULTS 
A. Field-induced AF long range order 

1. Order parameter and phase diagram 

The main manifestation of the field-induced transi- 
tion in IPA-CuCla is the onset of commensurate trans- 
verse long-range AF order. The high-field ordered 
phase is characterized by new Bragg reflections of type 
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FIG. 1: Field dependence of the (0.5, 0, 0) peak intensity mea- 
sured in at different temperatures (symbols). The lines are 
power-law fits to the data over a f T field range. Inset: The 
value of f3 fitted over a 1 T field range plotted as a function 
of temperature. The line is a guide for the eye. 



on each ladder rung are parallel to each other. The field 
dependencies of the (0.5,0,0) peak intensity measured 
at several temperatures are shown in the main panel of 
Fig. [T] (symbols). The order parameter critical index /3 is 
defined as ttiaf k [H -H^f, where mAF is the ordered 
AF staggered moment, and the magnetic Bragg intensity 
scales as / oc m 



AF- 



In Ref. [18| we demonstrated that 
power law fits over a progressively shrinking field range 
at T = 50 mK yield p close to 0.5. In the present work, 
power-law fits were performed over a fixed field range 
of 1 T above He at each temperature, and are plotted 
in solid lines in Fig. [1] The inset shows the thus mea- 
sured temperature evolution of (3. This critical exponent 
rapidly decreases upon heating, approaching the value 
/? — 0.25 at elevated temperatures. It appears that with 
both the fitting range and temperature approaching zero, 
/3 will indeed approach the Mean Field value of 0.5, as 
expected for a BEC-type transition. 

Another characteristic of the phase transition is the 
critical index v, defined as Tc oc {H — HcY . The mea- 
sured temperature dependence of the transition field is 
shown in Fig. [51 A power law fit to this experimental 
phase boundary gives v = 0.71(4). This value is close 
to = 2/3 expected for a 3D BEC of magnonSfi^ but 
should be regarded with some caution. Indeed, the den- 
sity of data points available at T — > is barely sufficient 
for an accurate estimate. 
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FIG. 2: Measured temperature dependence of the transition 
field in IPA-CuCls (symbols). The line is a power-law fit to 
the data up to T = 500 mK. 



{ "^^^^ ,k,l), with h, k and Z-integersji^. The magnetic 
structure at H — 12 T applied along the c direction 
was determined from 48 independent Bragg intensities 
measured using Setup VI. The AF ordered moment of 
0.49(1)/xb was found to be oriented almost parallel to the 
crystallographic b direction. As expected, the two spins 



B. Model dynamic structure factor 

Before discussing the results of inelastic measurements, 
we shall introduce an analytical model magnetic neutron 
scattering cross section for IPA-CuCla for use in the data 
analysis. A simple formula that accurately describes the 
dispersion of the de gene rate magnon triplet at 7? = 
was discussed in Ref. [l9l : 



{hw^°^f = {hLUofcos'^inh) + Alsm^{TTh) + 
+ CoSin^(27r/i). 



(1) 



The parameters were previously determined experimen- 
tally: hojo = 4.08(9) meV, Aq = 1.17(1) meV and 
Co = 2.15(9) meV. Equation [T] follows the structural pe- 
riodicity of the spin ladder: wave vectors with /i-integer 
and h- half-integer are not, in general, equivalent. For 
AF spin correlations h = 0.5 is the ID zone-center, while 
the zone boundary is located ai h — or h — I. Apart 
from the gap A, an important parameter is the velocity, 
that we defined as 



Vo = 



{0)^ 



A2 
I 



d{2TTh) 



1/1=0.5 



[(ftwo) 



A2]/4 = 2.9mcV. 



(2) 



In an external field the triplet spectrum is split due to 
Zeeman effect. For H < He the dispersion relation for 
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the i-th magnon branch with a field-dependent gap A^'^ 
can be written as: 



Ao. 



(3) 



Here we assume that the field only induces an overall en- 
ergy shift, without affecting the shape of the dispersion 
curve. Such behavior is indeed consistent with experi- 
ment (see below), and reflects the fact that the Zecman 
term commutes with the Heisenberg spin Hamiltonian. 

The dynamic neutron cross section for each mode was 
written as: 



dndE' 



(4) 



Here f{q) is the magnetic form factor for Cu^+ and Ai{q) 
is the overall intensity of mode i. The latter implicitly in- 
cludes the mode's structure factor S'i(q), the polarization 
factors for magnetic scattering of unpolarized neutrons, 
as well as the appropriate components of the g-tensor. 
The parameters of the model are the gap energies A.; 
and, at each given wave vector, the intensities Ai for 
each magnon branch. 

Equation 2] was also used to describe excitations in 
the high field phase. However, due to the onset of AF 
long-range order that doubles the structural period, at 
H > He h = and h = 0.5 become equivalent wave 
vectors, and the ID zone-boundary is shifted to h = 0.25 
or ft = 0.75. In this regime the dispersion relation for the 
i-th branch was written as: 



(ft^(' 



+ sin^(27r/i). 



(5) 



The parameters are the gap energies A^, velocities Vi = 
Ci, and, at any given wave vector, the intensities of each 
mode. 
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FIG. 3: Energy scan (raw data) collected in IPA-CuCls at 
H = 9 T using setup III, at q = (1.5,0,0) (solid symbols) 
and background measured away from the ID AF zone-center, 
at q = (1.3, 0, 0) and q = (1.7, 0, 0) (open symbols). 



the particular wave vector transfer. A similar procedure 
was applied to constant-g data at other wave vectors, or 
to those collected using most other experimental setups. 



Background-subtracted scans collected at using Setup 
III at various fields applied along the b axis are shown 
in Fig. dfa-d). As expected, the single peak seen at 
H = splits into three components at higher fields. 
The gap in the lower-energy mode extrapolates to zero at 
He ~ 9.5 T. Close to and above the critical field only two 
of the three peaks remain visible. In contrast with the 
previously studied anisotropic Haldane gap compound 
NDMAP, the gap in the lower-energy mode does not re- 
open above He- 



C. Spin hydrodynamics 

Inelastic data collected near the ID AF zone-centers 
defined by half-integer values of h provide information 
on the field dependence of the gap energies and the po- 
larizations of excitations. 



1. Gap energies 

A representative raw constant-q scan measured at the 
ID AF zone-center qi = (1.5,0,0) using Setup III at 
= 9 T is plotted in Fig. [3] in solid symbols. The 
background was measured at several fields as an average 
of scans at q'^ = (1.3,0,0) and q" — (1.7,0,0). Since no 
field dependence of the background was detected, it was 
further averaged to include all experimental fields (open 
symbols in Fig. [3]) and fit to a Gaussian centered at zero 
energy transfer plus a constant (solid line) . The resulting 
function was then subtracted from all scans measured at 



Very similar data at equivalent wave vectors were col- 
lected using Setup I for H\\a., and for H\\h using Setup II 
(Fig. O and Setup IV. In all cases, individual scans were 
analyzed using the model cross section described in the 
previous section. Eq. [5]was numerically convoluted with 
the 4-dimensional E — q resolution function of the 3-axis 
spectrometer, calculated in the Popovici approximation. 
At H > He, Eq. [5] was utilized whenever the entire dis- 
persion relation of each mode was measured, and thus 
both Ai and Ci could be independently determined (see 
below). In cases where the data collection was restricted 
to the ID AF zone-center, Eq. [S] was utilized instead, 
and only the gap energies and intensities were varied. 
For H < He Eq. [9] was used in all cases. Typical fits are 
shown in solid lines in Figs. S] and O In Fig. 2] the shaded 
areas represent partial contribution of each of the three 
members of the magnon triplet. The field dependencies 
of gap energies deduced from fits to the data collected us- 
ing all experimental conditions are summarized in Fig.[Sl 
In order to compensate for the g-tensor anisotropy, the 
abscissa in actually shows the reduced field value gH/2. 
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FIG. 5: Field dependence of background-subtracted inelastic 
scattering in IFA-CuCla measured at two almost orthogonal 
wave vectors q2 = (1.5,1.5,0) and qs — (1.5,-1.1,0) using 
Setup II. The reduced intensity at q2 is due to neutron ab- 
sorbtion in the sample. Solid lines are model cross section 
fits, as described in the text. 



FIG. 4: Field dependence of background-subtracted inelastic 
scattering in IPA-CuCla measured at the ID AF zone-center 
qi = (1.5,0,0) using Setup III. Solid lines are model cross 
section fits, as described in the text. Shaded areas are partial 
contributions of three separate excitation branches. 



2. Goldstone mode 



fits using Eq. [5] with A = 0, and assuming a constant 
background. Figure [7Jc) shows a cut through the TOF 
data of setup V taken between 0.8 meV and 1.2 meV en- 
ergy transfer. The sohd lines are empyrical Gaussian fits. 
From this analysis we obtain the spin wave velocities of 
the Goldstone mode vq = 1.74 meV and vq = 1.86 meV, 
for H = 11.5 T and = 14 T applied along the b axis, 
respectively. 



As mentioned above, aX H > the gap in the lower 
mode in IPA-CuCla does not re-open. However, this does 
not mean that the excitation itself disappears as, for in- 
stance, in the 5* = 1 dimerized chain system NTENP?^ 
In IPA-CuCla the lower mode persists above the critical 
field, but remains gapless, to within experimental resolu- 
tion. As discussed in Ref . [l^, this excitation becomes the 
Goldstone mode associated with the spontaneous break- 
ing of rotational symmetry due to magnetic ordering. It 
is best observed in g-scans at low energy transfers, as 
shown in Fig. [T] Here panels (a) and (b) display 3- 
axis data collected using Setup IV. The solid lines are 



3. Intensities and polarizations of excitations 

Since magnetic interactions along the crystallographic 
b axis are totally negligible, all wave vectors with the 
same values of h and I can be considered equivalent. In 
other words, the fc-dependence of magnetic scattering 
is fully decoupled from the energy dependence. Polar- 
ization information can be deduced from a comparison 
of constant-g scans collected at different yet equivalent 
wave vectors, thanks to neutrons being only scattered 
by fiuctuations of spin components perpendicular to the 
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gH/2 (T) 

FIG. 6: Gap energies measured in IPA-CuCls at T < 100 mK 
as a function of reduced field for different experimental ge- 
ometries (symbols). Solid and dash-dot lines are guides for 
the eye. 



momentum transfer. 

To take advantage of this selectivity, we compared 
scans collected at three equivalent ID zone-centers, sep- 
arated by roughly 45° in reciprocal space. In Setups I 
and III the data were taken at qi = (1.5,0,0) (Fig. [5]), 
to complement the scans measured with Setup II at 
q2 = (1.5,1.5,0) andqa = (0.5,-1.1,0) (Fig.[5l). These 
wave vectors form angles of ai = 0, as = 46°, and 
as = —45° with the crystallographic a* axis, respec- 
tively. Experimentally, for all geometries, the intensity of 
all gapped magnons remains constant as long as H < He, 
but decreases above the transition field. The intensity of 
the central magnon branch, deduced from model cross 
section fits to constant-q scans, is plotted against field in 
Fig. IHa). The data were normalized to the intensity of 
the single inelastic peak seen a.t H — 0. 

To within experimental statistics, the measured inten- 
sities of all visible gapped magnons remain proportionate 
at all three wave vectors in all fields. The corresponding 
intensity ratios are plotted in Fig. [Sl^b). In all cases the 
central mode remains roughly twice as strong the up- 
per and lower ones, and half as strong as the unresolved 
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FIG. 7: g-scans measured in IPA-CuCls at low energies in the 
BEC phase using Setup IV (a,b) and Setup V (c), showing 
the gapless Goldstone mode. In (a) and (b) the solid lines are 
fits to the data using a parameterized model cross section, 
as described in the text. In (c) The solid lines are empirical 
gaussian fits. 



triplet peak at H = 0. For H < He this qualitative 
picture is consistent with expectations. Indeed, the cen- 
tral mode is linearly polarized along the direction of ap- 
plied field (z-axis). It contributes equally to scattering 
at all wave vectors in the {x,y) scattering plane. The 
upper and lower modes are expected to have circular po- 
larization with Sz = ±1. The contribution of each one 
is evenly split between S'^'^ {q\\ , lu) and S'^^(g||,w). If x 
is chosen along the scattering vector, only the S"^"^ {q^^ , uj) 
contribution is detected, resulting in a two-fold reduction 
of intensity compared to linear polarization. A ratio of 
slightly larger than 1/2 is due to the anisotropy of the g- 
tensor in IPA-CuCla. The magnetic neutron cross section 
for a given spin polarization is proportional to the square 
of the corresponding g-factor, which accounts for an en- 
hancement of scattering in the {x, y) plane. A somewhat 
surprising experimental result is that the intensity ratio 
of the upper and central modes remains unchanged in the 
BEC phase. This behavior suggests that in experimen- 
tally accessible field, even at H > He where rotational 
symmetry around the z axis broken by long range order- 
ing, the intermixing of states with different polarization 
is not particularly large. 
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TABLE L Gap energies and spin wave velocities obtained 
from the measured magnon dispersion curves at three different 
fields. 
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FIG. 8: (a) Intensity of the central magnon branch measured 
in IPA-CuCla at T < 100 mK as a function of reduced field 
for different experimental geometries and scattering vectors 
(symbols). The intensity is normalized to that of the single 
peak measured at = 0. (b) Ratio of intensities of the 
upper and central (open symbols) and lower and central (dot- 
centered symbols) magnon branches measured in IPA-CuCls 
at T < 100 mK as a function of reduced field. In both panels 
the solid lines are guides for the eye. 



D. Spin dynamics away from the zone-center 

1. Magnon dispersion curves 

The dispersion of magnons across the Brilloin zone was 
measured in a series of constant-g scans. The data were 
taken at H ^ 9 T, H = 11.5 T and = 14 T using Se- 
tups IV and V. Due to geometric limitations on accessible 
energy transfers at different wave vectors, the dispersion 
relation of only the middle and lower modes were mea- 
sured in 3-axis experiments. These data were analyzed 
using the model cross section described previously. A 
semi-global fitting approach was implemented. Separate 
intensity prefactors were used for each individual scan, 
but the same dispersion relation was applied to all scans 
simultaneously at each field. The results of these fits is 
shown in solid lines in Fig. [9fa,c). Symbols represent 
fits to individual scans. In Fig. HJa) open symbols are 
data from Ref. [H collected at H = 0. For Fig. Mh), 
constant-q|| cutouts from the 2-dimensional TOF data 





H ^9T 


H = 11.5 T 


= 14 T 


Ai 


0.10(8) meV 


(fixed) 


(fixed) 


Aa 


1.17 meV (fixed) 


1.52(5) meV 


1.72(2) meV 


A3 


2.26(2) meV 


2.98(2) meV 




VI 


2.9 mcV (fixed) 


1.74(3) meV 


1.86(2) meV 


V2 


2.9 meV (fixed) 


2.29(6) meV 


1.84(3) meV 


V3 


2.9 meV (fixed) 


2.67(4) meV 





were fit to Gaussian profiles. The resulting dispersion 
relation (symbols) was then analyzed using Eq. O (solid 
lines) . The key parameters of the fits are summarized in 
Table U In all panels, a transition to dotted lines indi- 
cate that either the area is outside the experimental scan 
range or that no identifiable magnon peak was detected. 



2. Spectrum termination 

One of the main findings of Ref. was the discovery of 
an abrupt truncation of the magnon branch in IPA-CuCla 
at a critical wave vector h = ^ + he, he = 0.35. This phe- 
nomenon was attributed to 2-magnon decay, in analogy 
with 2-particle decay of rotons in liquid '*He. Similar 
physics was observed in the quasi-2D system PHCCi^ 
Further theoretical studie a^^'^^ showed that in one di- 
mension the effect of 2-particle decay is particularly dra- 
matic, and always leads to a total removal of the single- 
particle pole from the wave vector-energy domain of the 
2-particle continuum. One of the goals of the present 
work was to learn how the spectrum termination in IPA- 
CuCla is affected by an external magnetic field, particu- 
larly above He- 

In Fig. [To] we show a series of raw constant-g scans col- 
lected using Setup IV near the maximum of the measured 
dispersion, at = 9 T (just below the critical field). For 
/i < i -|- he, two peaks corresponding to the middle and 
lower magnon branches are clearly visible (Fig. llOb .b). 
The 3rd magnon is outside the scan range. The solid 
lines are fits using the model cross section given by Eqs.[4] 
and[9l convoluted with the spectrometer resolution. Both 
magnon peaks abruptly vanish at /i > ^ + he, as shown 
in Fig.fTOh. This behavior is quantified by the plot of the 
fitted intensity of the middle magnon branch, corrected 
for the magnetic form factor, against h in Fig. [Tl] (grey 
symbols). We conclude that at H = 9 T the spectrum 
termination is exactly the same as in zero field. Indeed, 
the trends visualized in Figure[Tn]are very similar to those 
shown for H = in Fig. 2(c,d) of Ref. ,19. 

The spectrum termination picture abruptly changes 
above the critical field. Figure [12] shows the h- 
dependence of constant-g scans collected at H — lA T> 
He using Setup IV. Only the "middle" branch of magnons 
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FIG. 9; Measured dispersion of magnons in IPA-CuCls below (a) and above (b,c) the critical field. Symbols are fits to individual 
constant-g scans, and lines are semi-global fits at each field, as described in the text. 
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FIG. 10: Constant-q scans measured in IPA-CuCla at H = 
9 T using Setup IV (symbols). Solid lines are fits to the 
data as described in the text. For the two visible branches of 
magnons the spectrum is terminated at ft = i + he, he ~ 0.8, 
as at H = 0^. 
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FIG. 11: Measured wave vector dependence of the intensity of 
the middle magnon branch at H = 9 T< He (grey symbols) 
and H — 14 T> He (solid symbols). Lines are guides for the 
eye. 



is visible, but does not extend beyond h^+h'^, h'^ = 0.25, 
which in the high field phase is the new ID magnetic 
zone boundary. The corresponding /i-dependence of in- 
tegrated magnon intensity is shown in Fig. [Tl] in solid 
circles. 

In Fig. [T^ . note that at ft, = 0.7 the observed peak is 
much narrower than those seen at = 9 T (Fig. [TOk). 
This, is entirely an effect of experimental wave vector 
resolution. All observed width are in fact instrumental 
at all fields. The narrow peak at = 14 T is due to that 
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fi(o (meV) 

FIG. 12: Constant-q scans measured in IPA-CuCls a.t H = 
14 T using Setup IV (symbols). Solid lines are fits to the 
data as described in the text. The magnon branch vanishes 
at ft = i + h'^,h'^ 0.25. 

h = 0.7 is close to the new zone-boundary, so that the 
local slope of the dispersion curve is shallow (Eq. [5]) . In 
contrast, at if < He, h — 0.75 is still far from where the 
dispersion levels out, and the local slope is steep (Eq.[9]). 
A steep slope leads to a broader peak, due to poor wave 
vector resolution, and a relatively good energy resolution. 
Similar arguments explain why at-ff = 9Tat/i = 0.8 
(Fig. 1 10b ) the magnons appear narrower than and h = 
0.7. 



3. Bandwidth collapse 

A very important result reported in Ref . [l^ is the dras- 
tic reduction of the bandwidth of the middle (and, prob- 
ably, top) magnon branch at high fields. From Fig.[S^ we 
see that there is practically no change in the dispersion of 
the middle magnon branch between H = Q and H — 9T, 
which means that the bandwidth collapse occurs only at 
H > He- The bandwidth continues to decrease in the 
entire accessible field range in the high-field state, from 
2.9 meV at ii = 9 T to as little at 0.8 meV at H = 14 T. 

A surprising finding is that between He and H = 
11.5 T the change of the spectrum of the central magnon 
is not continuous. In fact, there appears to be a field 
range in which large-bandwidth and low-bandwidth "ver- 
sions" of the middle excitation branch co-exist. This is 




ho) (meV) 

FIG. 13: A q = (0.7, -0.4, 0) scan measured in IPA-CuCls at 
H = 10.5 T using Setup IV (symbols). Three peaks are ob- 
served, and can be interpreted as a simultaneous observation 
of a 2-peak feature as in Fig. llOb (upper and lower peaks) 
and a single sharp peak as in Fig. 112b (central peak) . 



illustrated by Fig. [T3] that shows a scan taken at h = 0.7 
in an ii = 10.5 T external field. There clearly are three 
peaks in this scan. The two broad peaks can be associ- 
ated with the two highly dispersive magnons seen below 
He at H — 9 T (Fig. [TUb). At the same time, the sharp 
peak resembles that corresponding to the middle mode 
at higher fields (Fig. [HJi). 

The two versions of the middle magnon branch appar- 
ently coexist even at ii = 11.5 T, though the highly dis- 
persive component is considerably weakened. This can 
be inferred frorn_the false color plot of the TOF data in 
Fig. 1 of Ref. H^, where a portion of the middle mode 
seems to continue all the way to 4 meV as h 0.8, even 
as a stronger component levels off at 2.7 meV. Due to 
poor experimental resolution, low intensity, and a clear 
dominance of the narrow-bandwidth component, it is not 
possible to clearly identify the more dispersive compo- 
nent in constant-q cut-outs at this field. At = 14 T 
this excitation seems to have totally disappeared, and 
only the narrow-bandwidth component is detected. 



IV. DISCUSSION 

The behavior of spin correlations in IPA-CuCla under 
field can be summarized as follows. 

1) At -ff < He the system remains non-magnetic. 
All the way to the critical field the three lower-energy 
magnon branches are exact replicas of the degenerate 
magnon at H — 0, but are shifted in energy due to Zee- 
man effect, according to the spin projection that they 
carry: Sz = 0, ±1. This behavior is fully consistent with 
the Zeeman term commuting with the Heisenberg Hamil- 
tonian, and the ground state being an S" = singlet. 

2) To within experimental accuracy, the critical prop- 
erties of the transition at He and T — > are those of 
BEC in a dilute Bose gas. This is exactly what theory 
predicts in the case of ideal axially symmetryJ^i^i As 
the softening magnons condense in the ground state at 
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iJc, the initial rotational 0{2) symmetry of the system 
is spontaneously broken by the emerging AF long range 
order. This process is fully equivalent to the formation 
of a conventional Bose-Einstein condensate, where it is 
the C/(l) — 0(2) gauge symmetry that is being sponta- 
neously violated by the macroscopic phase of the conden- 
sate's wavefunction. In our case, at T = the density 
of magnons is strictly zero for H < He, and vanishingly 
small just above the transition. The phenomenon is then 
accurately described in the limit of negligible quasipar- 
ticle interactions. Deviations from this model become 
apparent with increasing temperature, due to a thermal 
excitation of magnons and a finite magnon density at the 
transition point. Experimentally, this thermal effect is 
particularly strong for the order-parameter critical index 
(3 that decreases rapidly with temperature. 

3) The BEC nature for the phase transition is consis- 
tent with the observed behavior of long-wavelength exci- 
tations near the AF point at H > He- The lower of the 
three magnons becomes the Goldstone mode, associated 
with the 0(2) symmetry breaking. To within experi- 
mental resolution it remains gapless at all fields above 
He- The intensity of this mode is decreased significantly 
compared to that below the transition field. At H > 
the gap in the middle magnon branch begins to increase 
linearly. Simultaneously, the rate of increase for the gap 
in the upper mode doubles. The intensities of both the 
middle and upper modes are suppressed with increasing 
field, while the velocities of all three excitation branches 
progressively decrease. 

This behavior is similar to that of other gapped quan- 
tum magnets undergoing field-induced ordering transi- 
tions, such as the S' = 1/2 dimer compound TlCuClar" 
the S = 3/2 material Cs3Cr2Brg^ or the Haldane spin 
chain system NDMAPj^i^ The key difference, of course, 
is that in those systems where magnetic anisotropy is 
present, the transition is of an Ising, rather than BEC 
universality class, and the lower mode re-acquires a gap 
in the ordered phase. ^ Field-theoretical descriptions of 
gapped spin chains^ii^^, as well as a bond-operator treat- 
ment of the coupled-dimer modelji^ attribute the changes 
of magnon gaps and velocities to an admixture of the 
higher-energy triplet modes to the ground state. This 
interpretation can be extended to our case of coupled 
spin ladders. In particular, at H > He rotational sym- 
metry is broken by the staggered moment, and Sz ceases 
to be a good quantum number for quasiparticles. The po- 
larizations of all magnon branches must become mixed. 
The Goldstone mode is thus distinct from the Sz = I 
magnon at H < He, just as the two gap modes above 
the condensate can no longer be identified with Sz — 
and Sz — —I quasiparticles at low fields. For the acces- 
sible range of magnetic fields, this polarization mixing 
in IPA-CuCls can not be too strong though, as it was 
not directly detected in our polarization studies. This is 
undoubtedly due to the fact that even at 7J = 14 T the 
ordered moment is still far from its saturation value. On 
the other hand, the method used is, admittedly, not the 



most accurate one: small shifts in polarization will be 
easier to detect in future polarized neutron experiments. 

4) Short wave length spin dynamics of IPA-CuCla at 
H > He IS rather complex. The Goldstone mode is not 
visible at higher energy transfers, away from the zone- 
center, presumably due to its low intensity. The middle 
magnon branch is actually replaced by two excitations 
that merge at the zone-center, but are distinct near the 
zone-boundary. One of these retains the steep dispersion 
relation seen at H < He- This component progressively 
weakens with increasing field and practically vanishes at 
_ff « 12 T. Simultaneously, the 2nd component gains in 
intensity. It has a much shallower dispersion, and entirely 
dominates at high fields. The bandwidth of this compo- 
nent progressively decreases with increasing H: the gap 
at the zone-center increase, and zone-boundary energy 
goes down. The limited data that has been accumulated 
on the 3rd (highest-energy) magnon branch suggest that 
it's behavior is similar to that of the middle mode. 

In our previous brief report the drastic restructuring 
of the excitation spectrum above He was attributed to a 
spontaneous breaking of a discrete microscopic symmetry 
at the transition point. — Unlike in a conventional BEC 
transition that breaks a continuous U{1) = 0(2) sym- 
metry, the transition in IFA-CuCla in addition to that 
doubles the magnetic structural period. The result is a 
folding of the Brillouin zone and an opening of anticross- 
ing gaps at the new zone boundaries {h — 0.25, h — 0.75, 
and equivalent positions), where each branch interacts 
with its own replica from an adjacent zone. This, in turn, 
leads to an abrupt reduction of the zone-boundary energy 
for the visible (lower) segments of the gapped magnons. 
The additional breaking of translational symmetry is not 
an inherent property of all magnon condensation transi- 
tion, but specific to IFA-CuCla and select other systems. 
For example, in the dimer material TlCuCls the peri- 
odicity of the magnetically ordered state is the same as 
that of the underlying crystal lattice. Correspondingly, 
the changes in the excitation spectrum at the transition 
point^^ arc continuous and much less drastic than in IPA- 
CuClg. 

The doubling of the ground state's translational pe- 
riodicity is also the primary cause for the shift of the 
critical wave vector of spectrum termination at H > He- 
As discussed in detail in Ref. 19,, at H = the magnon 
branch is truncated at a point beyond which two-particle 
decay becomes allowed by energy and momentum conser- 
vation. The particular shape of the magnon dispersion 
curve in IPA-CuCla is such that the critical wave vector 
is he ~ 0.35. Even though the magnon energies shift at 
< H < He, the termination point for each branch re- 
mains unchanged in that regime. Due to conservation 
of spin projection, the 5*^ = magnon, for example, 
can only decay into quasiparticles with 5*2 = —1 and 
5z = 1. Since both the initial and final states are eigen- 
states of the Zeeman term, none of the matrix elements 
are affected by the applied field. Moreover, the Zeeman 
shifts for the decay products are opposite and fully cancel 
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each other, so any energy and momentum conservation 
requirements remain as at if = 0. Similarly, below 
the field has no effect on the decay of the ^2 = ±1 par- 
ticles. 

This picture has to change qualitatively at ff > He- 
The breaking of rotational symmetry in the ground state 
removes the Sz conservation requirement. It is easy to see 
that in this case 2-particlc decay becomes possible at any 
wave vector beyond the new zone-boundary h = 0.75, so 
that h'^ = 0.25. The final state is a combination of the 
same type of magnon with a smaller wave vector, and a 
Goldstonc mode. 

Certain important details regarding the short wave- 
length spin dynamics of IPA-CuCls above the critical 
field are currently not understood. In particular, we have 
no obvious explanation for the apparent co-existence of 
wide- and narrow-bandwidth magnons just above He- 
Further theoretical studies of the problem are clearly in 
order. 



V. CONCLUSION 

To date, most studies of field-induced ordering tran- 
sitions in quantum magnets were focused on long wave 



length spin correlations. Indeed, these properties are di- 
rectly relevant to the phase transition itself. However, 
it appears that a lot is to be learned and understood 
about short wavelength spin correlations as well. This 
is particularly true when field-induced BEC of magnons 
is accompanied by a breaking of discrete translational 
symmetry, as in IPA-CuCls. 
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